Localized 4/ States and Dynamic Jahn- Teller Effect in PrC>2 
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Neutron spectroscopic measurements of the magnetic excitations in PrC>2 reveal (1) sharp peaks 
characteristic of transitions between levels of the 4/ 1 configuration of Pr 4+ split by the cubic crystal 
field, and (2) broad bands of scattering centred near 30 meV and 160 meV. We present a simple model 
based on a vibronic Hamiltonian that accounts for the main features of the data. The analysis shows 
that 90 ± 10% of the Pr ions have a localized 4/ 1 configuration and provides strong evidence for a 
dynamic Jahn- Teller effect in the V$ electronic ground state. 

PACS numbers: 71.28.+d, 71.70.Ej, 71.70.Ch, 78.70.Nx 



Anomalous static and dynamic magnetic properties of 
rare earth and actinide compounds containing strongly 
correlated electrons are usually attributed to Kondo or 
valence fluctuations, or to the formation of heavy fermion 
states. If there are unquenched orbital degrees of free- 
dom, however, then other mechanisms for anomalous be- 
havior are possible. One such mechanism is the coupling 
of electronic charge fluctuations to phonons, and leads 
to magnetoelastic and Jahn- Teller phenomena. In the 
case of localized electrons the dynamic Jahn- Teller ef- 
fect (DJTE) may quench electronic degrees of freedom, 
causing a reduction in the magnetic moment and forming 
states of mixed magnetic and phonon character |jj . 

Electron-lattice effects of this kind are of course well 
known but the absence of a clean reference system has 
undoubtedly hindered their identification among other 
anomalous properties, especially in Ce and U compounds. 
In this Letter we show that PrC>2 represents such a model 
DJTE system. We find that magnetoelastic coupling has 
a dramatic influence on the static and dynamic properties 
of PrC>2, and we provide the first direct measurements of 
the magnetovibrational spectrum in a rare earth DJTE 
system. We also show that the Pr ions exist almost en- 
tirely in a localized 4/ 1 configuration, in disagreement 
with an intermediate valence model for PrC>2. 

PrC>2 is an insulator with the fluorite structure type. 
It exhibits type-I antiferromagnetic ordering below Tn = 
14 K with an anomalously low ordered moment of 
^0.6 /ib Hi- The electronic ground state of PrC>2 has, 
together with CeC>2 and TbC>2, been controversial for 
many years Core-level photoabsorption and pho- 

tocmission spectra have been interpreted by some au- 
thors |j]]4|] in terms of an intermediate valence ground 
state comprised of a roughly 50:50 admixture of local- 
ized (i.e. atomic- like) 4/ 1 and 4/ 2 states. Others H|), 
however, believe in a localized 4/ 1 configuration with a 
degree of covalent mixing such that the oxygen 2p va- 
lence band contains some extended states of / symme- 



try. Publication of an influential model (jj describing 
an intermediate valence state in the anomalous cuprate 
PrBa2Cu306+:r (PrBCO) similar to that proposed for 
the lanthanide dioxides has increased interest in Pr02, 
making it a key reference compound in high energy spec- 
troscopic studies of PrBCO || and related materials. 

The aim of the present work was to clarify the elec- 
tronic ground state of Pr02- Neutron spectroscopy 
probes electronic transitions within the ground state con- 
figuration both within and between the 2S+1 Lj terms, 
including splittings due to the crystalline electric field 
(CEF) or exchange fields. The low neutron energies 
(< 1 eV) and weak coupling to the electronic angular 
momentum ensure that final state effects are insignifi- 
cant. The experiments were performed on the High En- 
ergy Transfer (HET) chopper spectrometer at the ISIS 
Facility. On HET neutrons of well-defined incident en- 
ergy Ei are delivered to the sample in short pulses. Spec- 
tra are recorded as a function of flight time in banks of 
detectors surrounding the incident beam direction. Be- 
cause the spectra are recorded at constant scattering an- 
gle (j> the scattering vector Q varies with energy E thus: 
h 2 Q 2 /2m = 2Ei-E - 2cos<j)y/E i (E i - E). 

Polycrystallinc samples of Pr02 were prepared by ox- 
idation of commercially obtained PrgOn. The starting 
material was baked in air at 1000 °C for several hours 
and then annealed either in flowing O2 at 280 °C for ap- 
proximately 20 days or under 300-500 atm pressure of O2 
at 350 °C for 3-5 days. The products were checked by x- 
ray diffraction and no trace of residual P^On could be 
detected. From this, and a comparison of our neutron 
data with that of PrgOn we estimate an upper limit 
of 1% on the amount of PrgOn impurity in our samples. 

For the neutron measurements approximately 10 g of 
Pr02 was wrapped in aluminium foil and mounted in con- 
tact with the cold head of a closed-cyle refrigerator. Data 
were collected with E\ = 30, 180, 450, 750 and 1200 meV, 
and the scattering from vanadium was used to normalize 
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the spectra and to convert them into units of absolute 
scattering cross-section. The presented spectra corre- 
spond to the partial differential cross-section d 2 a/dfldEf 
multiplied by the factor fcj/fcf ]To| | , where fci and kf are 
the initial and final neutron wavevectors and Ef is the 
final energy. Measurements were made under identical 
conditions on a 10 g sample of polycrystalline CeC>2 in 
order to estimate the non-magnetic (phonon + multiple 
scattering) background in the PrC-2 data. CeC-2 is an 
almost ideal material for this purpose as it has the same 
structure as PrC>2 and the scattering cross-sections of Pr 
and Ce differ by only 10%. The estimated transmission of 
the PrC>2 and CeC>2 samples was 93-97% (depending on 
E and E{), and the spectra were corrected accordingly. 

Fig. 1(a) shows spectra of PrC>2 and CeC>2 with E\ = 
180 meV averaged over the angular ranges 9-19° ((</>) — 
14°) and 125-139° ((0) = 133°). At high angles the 
scattering is dominated by one-phonon processes, and 
so the 133° data measures the phonon density of states 
weighted by the neutron cross-section of the atoms. The 
133° data for the two compounds arc similar except that 
some PrC>2 peaks are shifted to lower energies relative to 
CeC>2, most noticeably below 30 meV. At low angles more 
marked differences are observed. The PrC>2 data contains 
a sharp peak at 131 meV and a broad peak centred on 
~30meV extending from lOmeV to 80meV, whereas the 
CeC>2 spectrum shows little structure and lies below the 
PrC>2 spectrum at all energies. The 131 meV peak has 
been observed previously || . Both the broad and sharp 
peaks in the PrC>2 data decrease in intensity with increas- 
ing Q consistent with a typical 4/ magnetic form factor, 
and so we identify these features as magnetic in origin. 
Conversely, the Ce02 signal systematically increases in 
intensity with Q, and this implies that the Ce02 scatter- 
ing is non-magnetic as expected for a 4/° configuration 
with no admixture of localized 4/ 1 states. 

Several other magnetic features were observed in the 
PrC>2 data in addition to those already described. The 
runs with E[ = 750 meV and 1200 meV detected two 
more peaks. These are shown in Figs. 1(b) and (c). The 
peak centred near 350 meV in the 750 meV run is sig- 
nificantly broader than the resolution, whereas the peak 
centred at 730 meV in the 1200 meV run is resolution- 
limited. The 450 meV run revealed a shoulder of scatter- 
ing above the 131 meV peak, centred at ~160meV. This 
feature is displayed in the inset to Fig. 1(d). Finally, 
the Ei = 30meV run revealed a peak centred at 3meV 
shown in Fig. 1(d). This peak shifted to lower energies as 
the temperature was raised, becoming quasielastic above 
~15 K. Given that Tn = 14 K we can confidently at- 
tribute this 3meV peak to spin wave excitations of the 
antifcrromagnetically ordered ground state. We found 
no evidence in our data for an excitation at 5.4 meV re- 
ported in a previous work [ fT2[ , and so we conclude that 
this feature is not intrinsic to Pr02- 

Fig. 1(d) shows data from several runs after correction 



for the non-magnetic scattering and for the calculated 
(see below) magnetic form factor. The form factor cor- 
rection extrapolates the data to zero Q and hence allows 
comparison of runs with different incident energies. 

We begin our interpretation of the results by consider- 
ing the effect of the CEF. The splitting of a localized 4/ 1 
configuration in a cubic CEF is described by the spin- 
orbit coupling constant £, which sets the 2 F 5 / 2 ~ > 2 F-j/ 2 
separation, and the fourth and sixth order CEF param- 
eters Bq and Bq. Using the free ion value of £ for Pr 4+ 
JLlj Kern et al |0] have calculated the energy levels as a 
function of the ratio Bq/Bq. They assigned the 131 meV 
peak to the T% — > T-? transition of 2 Fc l / 2 , and by esti- 
mating the ratio Bq /Bq from a point charge model they 
arrived at a parameterization of the CEF. Kern et a/'s 
model predicts the levels in the CEF-split 2 F 7 / 2 term to 
be ~320meV (r' 6 ), -390meV (1^) and ~580meV (r' 7 ) 
above the ground state in rough accord with the peaks 
in Figs. 1(b) and (c). We conclude, therefore, that these 
high energy peaks arise from the 2 F 5 / 2 — ► 2 F 7 / 2 intermul- 
tiplet transitions, the only proviso being that the signal in 
Fig. 1(b) must encompass both the T' 6 and T' 8 levels. As- 
suming this to be the case we refined the CEF model with 
the 14 states | J; mj) of 2 F^/ 2 and 2 F 7 / 2 as a basis. With 
all 3 parameters allowed to vary independently the fit- 
ting procedure converged yielding Bq — (—776 ± 8) meV, 
Bl = (207 ± 5) meV and £ = (100.5 ± 1) meV [jj and 
gave excellent agreement with the observed energies. 

With the eigenfunctions from the CEF refinement we 
can now evaluate the cross-sections for the observed tran- 
sitions as a function of Q. The formulae [jl0| contain 
radial moments of the 4/ wavefunction, and values for 
Pr 4+ were taken from Ref. pl|. The curve in Fig. 1(b) 
represents the calculated cross-sections for the Tg — > T' 6 
and Tg, — > T' 8 transitions, and the curve on Fig. 1(c) is 
likewise for the Tg — > T' 7 transition. Considering that no 
additional parameters are needed to obtain the calculated 
cross-sections the level of agreement with the data is very 
good. The same cannot be said, however, when the cross- 
sections for the Tg — > Tg and Tg — > T 7 transitions are cal- 
culated. At zero Q these are 182 and 99mbsr~ 1 Pr~ 1 re- 
spectively fl6|l , considerably greater than the integrated 
spectral weights 47 and 48 mb sr~ 1 Pr _1 in the 3 meV and 
131 meV peaks shown in Fig. 1(d) |T^| . 

We will now describe a model that both accounts for 
the puzzling intensities just mentioned and explains the 
origin of the 10-80 meV broad peak and the 160 meV 
shoulder. Our hypothesis is that these effects derive from 
a strong coupling between 4/ 1 electronic states and local 
dynamic lattice distortions. This coupling mixes elec- 
tronic and phonon degrees of freedom and causes a dy- 
namic Jahn- Teller effect in the Tg ground state. 

The existence of a DJTE in Pr02 was in fact sug- 
gested some time ago to explain the factor ^2 reduction 
in ordered moment relative to a pure Tg ground state 
H . An analogous moment reduction in UO2 had already 
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been attributed to a DJTE [Q, and so its occurrence 
in PrC>2 would not be a surprise. Moreover, the neutron 
data in Fig. 1(a) showing substantial phonon shifts in 
PrC>2 relative to CeC>2 provide direct evidence for a sig- 
nificant magnetoelastic interaction. To find out how this 
interaction influences the magnetic excitation spectrum 
we describe a simple model for the coupling of a set of 
CEF levels with phonons. Group theory shows that the 
T 8 CEF ground state couples in first order to local dis- 
tortions of either T3 or T§ symmetry, and point charge 
calculations indicate that the coupling strengths are com- 
parable H). For the sake of simplicity, however, we will 
restrict our model to three degenerate local distortions 
of T§ symmetry with vibrational frequency uj p h. We de- 
scribe the magnetoelastic coupling with the Hamiltonian 

H M E = g^ia, + a\)O u (1) 

i 

where g is the coupling constant, a, and a\ are phonon 
annihilation and creation operators, and Oi are the T§ 
symmetry quadrupolar operators: 0\ — (J y J z + J z J y )/2, 
2 = {J z Jx + J x Jz)/2, and O3 = (J x Jy + J y J x )/2. To in- 
vestigate the ordered magnetic moment we also included 
a molecular exchange field interaction H^ x = H • J. The 
magnitude of H was chosen so as to reproduce the ob- 
served 3 meV splitting of the ground state. 

We take as a basis the 24 states \4> n ), n = 1, . . . , 24, 
represented by |r 8 ;0), T 7 ; 0), |r 8 ; 1«), and |r 7 ;lW). 
These are products of the CEF eigenstates and the three 
T$ phonon modes (i = 1, 2, 3) containing either or 
1 quanta hujph- The energy separation A of iTgjO) and 
|T7; 0} is chosen so as to reproduce the observed peak 
at 131 meV. The vibronic states are obtained by diago- 
nalisation of the matrix (</>„/ \Hme + H-E X \(f> n ), and the 
magnetic part of the neutron cross-section is calculated 
from the matrix elements of J. 

We evaluated the model for a range of parameters. 
Small values of huiph are effective at mixing |rg;0) and 
rg;!*- 1 ') states creating a DJTE. The consequences are 
a transfer of intensity from the ground state into the 
phonon-like states and a reduction in the ordered mo- 
ment. A large TiuJ v h tends to create vibronic states in- 
volving the r 7 CEF excitation. To illustrate these ef- 
fects we plot on Fig. 2 the zero-Q energy spectrum for 
two parameter sets, one (I) with small and the other (II) 
large fiujp^. In both cases new peaks appear between the 
3meV and 131 meV levels. With set I there is a factor 
0.86 reduction in the ordered moment (independent of 
the direction of H) and a similar decrease in the inten- 
sity of the 3 meV peak. Set II produces a smaller effect on 
the ordered moment and 3 meV peak, but a new feature 
is the signal above the 131 meV peak. 

These results give us confidence that magnetoelastic 
coupling is responsible for the unusual energy spectrum 
and reduced moment of Pr02. Because of its simplicity 



we cannot expect a single-frequency model to match the 
data perfectly. In reality there will be coupling to local 
dynamic distortions with T3 as well as T§ symmetry, and 
these distortions will exist over a range of frequencies due 
to dispersion. The extent of the observed broad scat- 
tering is indicative that many frequencies are actually 
involved. A realistic model would also need to include 
states with >1 phonons. Nevertheless, the qualitative 
agreement achieved with the present model is satisfying. 

Before concluding we will address the nature of the / 
states. Our analysis has assumed a localized 4/ 1 configu- 
ration. We can check this assumption two ways. Firstly, 
we can look for the signature of localized 4/ 2 states in 
the range 200-300 meV where the 3 #4 — > 3 H$ transition 
occurs Q. No peaks are observed. Second, we can di- 
rectly determine the number of 4/ 1 states per Pr from 
the neutron cross-section via the sum rule for transitions 
within a J multiplet |Tcfl . The measured zero-Q cross- 
section integrated up to 220 meV, including the elastic 
scattering fll7|| , is 240mbsr~ 1 Pr~ 1 . This compares with 
the calculated value of 182 + 99 = 281mbsr" 1 Pr" 1 for 
the J = I multiplet Thus, by including the broad 
scattering in the integral we recover almost all the miss- 
ing intensity. Allowing a 10% uncertainty in the absolute 
calibration we conclude that 90 ± 10% of the Pr ions have 
a localized 4/ 1 configuration. The intermultiplet transi- 
tions shown in Figs. 1(b) and (c) provide further support 
for this conclusion since the calculated cross-sections for 
4/ 1 are in good agreement with the measured intensities. 

In summary, our findings support a tetravalent model 
for Pr02 with an almost fully occupied atomic 4/ orbital. 
The neutron data allow the possibility of ^10% under- 
occupancy relative to a pure 4/ 1 configuration, and this 
may indicate a degree of covalency. Any other occupied 
states of / symmetry must then be extended states in the 
O 2p valence band. These results rule out the strongly 
intermediate valence model for Pr02, and we anticipate 
that a similar picture will apply to Ce02 and TbC>2. We 
find compelling evidence that magnetoelastic coupling in 
PrC>2 creates a ground state with mixed electronic and 
vibrational degrees of freedom by means of a DJTE. Be- 
ing an insulator with a simple structure, localized / elec- 
trons and a large CEF splitting, PrC>2 would seem an 
ideal model system for further studies of the DJTE. 
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FIG. 1. Neutron inelastic scattering from Pr02 and Ce02 
measured at a temperature of 10 K. In (a) the incident neutron 
energy Ei was 180 meV, and spectra recorded at low and high 
scattering angles are shown, (b) and (c) show 2 F 5 / 2 — > 2 F 7 / 2 
intermultiplet transitions measured with Ei = 750 meV and 
1200 meV respectively. The solid lines are cross-sections cal- 
culated from the CEF model described in the text with widths 
corresponding to the instrumental resolution, (d) shows data 
corrected for the non-magnetic background and the Q depen- 
dence of the magnetic cross-section as described in the text. 
Main frame: E\ = 30 and 180 meV. Inset: Ei = 450 meV 
showing the shoulder to the 131 meV peak. 

FIG. 2. Zero-Q cross-section of Pr02 calculated from the 
magnetoelastic model described in the text with parameter 
sets I (TiWph = 12.5 meV, g = 8meV and A = 109 meV) and 
II (huj ph = 42meV, g = 9meV and A = 101 meV). In both 
cases H was 0.5 meV || (1, 1, 1). The peak widths correspond 
roughly with the experimental resolution in Fig. 1(d). 
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